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Abstract — Cooperation between terminals has been proposed 
to improve the reliability and throughput of wireless commu- 
nication. While recent work has shown that relay cooperation 
provides increased diversity, increased multiplexing gain over 
that offered by direct link has largely been unexplored. In this 
work we show that cooperative multiplexing gain can be achieved 
by using a half duplex relay. We capture relative distances 
between terminals in the high SNR diversity multiplexing tradeoff 
(DMT) framework. The DMT performance is then characterized 
for a network having a single antenna half-duplex relay between a 
single-antenna source and two-antenna destination. Our results 
show that the achievable multiplexing gain using cooperation 
can be greater than that of the direct link and is a function of 
the relative distance between source and relay compared to the 
destination. Moreover, for multiplexing gains less than 1, a simple 
scheme of the relay listening 1/3 of the time and transmitting 2/3 
of the time can achieve the 2 by 2 MIMO DMT. 

I. Introduction 

There is a growing interest in the design of cooperative 
schemes that provide diversity and multiplexing gain for com- 
munication via wireless relays. Cooperative diversity refers to 
the additional diversity gain (compared to direct link) offered 
by cooperation. Similarly if a relay provides additional degrees 

■ of freedom (compared to direct link) it is said to provide a 

\ cooperative multiplexing gain 

Diversity multiplexing tradeoff (DMT) |2| has been widely 
used to analyze and compare the performance of coopera- 
tive schemes. DMT for the half-duplex single relay network 
has been studied extensively in literature |3||4||5|. For the 
case with single antennas at all terminals the 2 x 1 MISO 
DMT bound has recently been shown |5| to be achievable. 
In this paper we study the DMT for the multiple antenna 
half-duplex relay channel having m, n and k antennas at 
source, destination and relay respectively. This was studied 
in IJl but results were shown only for the special case 
m = n = 1. We calculate the maximum achievable DMT 
for the m = l,n = 2, /c = 1 configuration. In the process 
we also demonstrate techniques that enable results for general 
m, n and k. m = l,n = 2,/c = 1 is the simplest configuration 
where relay cooperation provides additional multiplexing gain 
compared to direct link. We show that if source and relay are 
relatively close to each other, cooperative multiplexing gain is 
achievable even with half duplex relaying. Note that the full- 
duplex case has been studied in 1 1 1. Moreover, for multiplexing 
gains less than 1, if source-relay SNR (measured in dB) is at 
least two times the source-destination SNR a simple scheme 




Fig. 1. Relay channel with 2 antennas at destination and S ^ R proximity 
gain 7]. 



with the relay listening | of the time and transmitting | of 
the time can achieve the 2 x 2 MIMO DMT. 

These results lend fresh insight into the fundamental limits 
of cooperative multiplexing in the half-duplex relay channel. 
We demonstrate the use of two key techniques that enable our 
results. 

1) Distance between terminals: In most results it is seen 
that relative distances between source, relay and des- 
tination do not affect DMT performance of the relay 
channel. Since DMT is calculated at high SNR the path 
loss and therefore distances are not easily captured in 
results. We overcome this apparent limitation by scaling 
the average SNR's of the various links differently. 
Our approach enriches the DMT framework by adding 
insights about network geometry. 

2) MIMO with half duplex antenna: The min-cut capacity 
bound has been used in O to calculate an upper 
bound for DMT performance. Notice in Fig. [T] that 
the {S,R},{D} cut corresponds to a 2 x 2 MIMO 
system with one source antenna that remains active only 
for a fraction of total communication time (R is half 
duplex). It was noted |4| that an upper bound for mutual 
information across such a cut is hard to compute. Due 
to this, DMT bounds have only been reported for the 
special case of m = n = 1. 

In Sec IIII-AI we demonstrate a simple channel decom- 
position that allows us to compute the cut-set DMT 
bound for the m = l,n = 2, /c = 1 configuration. 
The technique can be applied towards computing DMT 
bounds for general m, n and k. Recent results in |7J O 
show that a simple relaying scheme called "quantize- 
map" can achieve a rate within constant gap of the cut- 
set capacity. In SecjlVlwe discuss this scheme and show 
that it achieves the cut-set DMT bound. 
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Fig. 2. d(r) for various values of proximity gain 77. 



II. System Model 



Consider the system in Fig. [T] with source S, relay R and 
destination D having 1 , 1 and 2 antennas respectively. Let 
Dj^j G {1, 2} denote the jth antenna at D. The channel gain 
forS^R is hsr, gains for S Dj are hgj and R Dj are 
hdj . All the channel gains are assumed to be flat fading having 
i.i.d. CA/'(0, 1) distribution. We assume quasi-static fading, 
i.e. once realized, channel gains remain unchanged for the 
duration of the codeword and change independently between 
codewords. Noise at all receivers is additive i.i.d. CA/'(0, 1) 
and independent of all other variables in the system. Transmit 
power at S and R is limited by an average power constraint. 
Since noise power at receiver is normalized to 1, the transmit 
power constraint is specified by the average Signal to Noise 
Ratio (SNR). R is assumed to operate under a half-duplex 
constraint. For simplicity it is assumed that transmissions at 
S and R are synchronous at symbol level. 

We assume an asymmetrical network geometry. S and R 
are modeled to be close to each other as compared to {S^ R} 
and D. S ^ D and R ^ D diXQ assumed to have the same 
average SNR denoted by p. S ^ R on the other hand is 
modeled to have SNR higher than p by a factor r] on dB 
scale, i.e. S ^ R average SNR is . The S ^ R channel 
(cooperation link) thus has r] — 1 more degrees of freedom 
than other channels in the network. We call 77 the proximity 
gain and assume r]>l. 

No channel state information (CSI) is available at S i.e. 
only average channel statistics p, r] are known. However, at D 
all channel realizations hsr, hsj^ hdj are completely known. 

We identify three models for relaying strategy. 

• Global: The relay uses knowledge of all instantaneous 
channel realizations to optimize its strategy. 

• Local: The relay can measure hsr and uses only this 
(local) information. 

• Blind: The relay only uses average channel statistics. 

The global strategy is discussed in Sec Jill while local and 
blind are discussed in Sec (V] 



III. Diversity-Multiplexing Tradeoff 

Theorem 3.1: The maximum achievable DMT for network 
described in Sec HIl is given by, 

{m\n{r] + 2, 4} - 3r < r < 1, 7^ > 1 

(27/ - 7/r - l)/{r] -r) l<r<2-^,7/>2 
^-2^ l<r<2-i l<7/<2 

(1) 

Corollary 3.2: For system model described in Sec Ull the 
maximum achievable multiplexing gain r* = mfr>o{r\d{r) = 
0} is, 

r* = 2 - - (2) 
V 

For a symmetrical geometry with all channels having the 
same degrees of freedom (r] — 1) we get r* = 1 i.e. 
cooperation doesn't provide additional maximum multiplexing 
gain. To enable higher multiplexing gain the 5 ^ channel 
(cooperation link) needs to have more degrees of freedom than 
the S ^ D channel (communication link). 

Let d2x2{r) represent the DMT of the 2 x 2 MIMO channel. 
For finite r] it can be seen that d{r) < d2x2{r) with strict 
inequality over a non-empty region of r. This suggests that 
for distributed antennas the finite capacity of the cooperation 
channel {S ^ R) poses a fundamental limitation on the 
achievable DMT performance. It can easily be verified that, 

lim d{r) d2x2{r) 
77^00 

Fig. [2] shows d{r) for several values of r]. 

We prove Theorem 13. II in two steps. In Sec IIII-Al we show 
that the cut- set DMT upper bound for network in Sec [III is 
given by ([T]). In Sec|IV]we show that this bound is achievable. 

A. Cut-Set DMT upper bound 

Let /(O < / < 1) denote a listen-transmit schedule for the 
half duplex relay. R listens for a fraction / (listening phase) 
of total communication time and transmits for fraction (1 — 
/) (cooperation phase). The two cuts of the network Cd = 
{S,R},{D} and Cs = {S}{R,D} are shown in Fig. [3] for 
these two phases. In the listening phase let denote the 
sequence of symbols transmitted by S while Y^^ and 
denote received signals at R and D respectively. Similarly for 
the cooperation phase X| and X^^ are the symbol sequences 
transmitted from S and R while Yj^ is received at D. The 
instantaneous mutual information across the two cuts can be 
written as, 

Ics = //(X^;YH,Yi,|X^) + (l-/)/(X|;Y?,|XH)(3) 
Icn = //(X^;Yi,) + (l-/)/(X|,XH;Y?,) (4) 

To maximize these mutual information expressions we 
need to choose zero-mean complex Gaussian distributions 
for X^,X| and X^^ that have covariance matrices which 
satisfy their respective average power constraints. Using these 
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Fig. 3. Two cuts of network during listen (/) and cooperation phase.(l — /). 
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Fig. 4. Parallel Channel Model for 2 x 2 MIMO. 



interference cancellation receiver achieves the mutual infor- 
mation of the MIMO channel. For this scheme gf and ^2 can 
be calculated to be, 



9l 



9l 



where h^^^ and h^yg respectively denote the perpendicular 
and parallel components of with respect to h^. 

Note that while g\ and ^2 correlated, hg, h^^g and h^y^ 
are mutually independent. The correlation between gl and gl 
can therefore be explicitly calculated. The destination decodes 
Xi^ in the presence of interference from X|. It then cancels 
Xi^ from its received signal before decoding X|. Therefore 
S effectively sees an interference free channel (with gain gi) 
to D during both listen and cooperation phases. 

I'c, = flog{l+pgf) 

+ (1 - /)[log(l + pgl) + log(l + pgj)] 
= \ogil + pgl) + il-f)\ogil+pgi) 



distributions we can write mutual information upper bounds 
for Ics and /^^ respectively. 



where = 



f\0g(l^p^\hsr?^p\\K\?) 

+ (l-/)log(l + p||h,|p) 

/max{log(l + p^|/i,,P),log(l + p||h,|p)} 

+ (l-/)log(l + p||h,|p) 

/log(l + p||h,|p) 

+ (1-/) logdet(I + pHHt) 



hsi 

hs2 



hrl 
hr2 



and H = [ hg ]. 



It can be verified that the approximation is tight within one 
bit. 

Note that the expression for /^^ is a linear combination 
of the capacities of two Raleigh fading Gaussian channels 
having correlated channel matrices Hi = [ hg ] and 
H = [ hg ]. Outage probability for the MIMO channel 
was calculated in O using eigenvalue decomposition of the 
channel matrix. Following the same technique there would 
require computing the joint eigenvalue distributions for the 
two correlated hermitian matrices, HiHJ and HH^ It was 
noted in f4\ that this is hard to compute. We propose an 
easier decomposition to solve this problem. The second term 
in I^^ is the capacity of a 2 x 2 MIMO channel which can 
be represented by two parallel Gaussian channels having gains 
gi and ^2 shown in Fig. (H The channel can be written as. 
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where = 1, wi^W2 ~ CA/'(0, 1). The capacity for 

Xi Vi is given by log(l + P9i)- It was shown in ||6l||2l that 
a D-BLAST transmission scheme with a MMSE successive 



Let asr,<^i and 0^2 represent channel realizations via the 
following variable transformations. 

log(l + p^|/i,,p) 



a. 



ai 



OL2 



lim 



logp 
Ijj^ log(l + /^fli) 

p^oo log p 

^.^ iog(i + pgi) 

P^oo log p 



This gives us simplified expressions for mutual information 
upper bounds. 



logp 



^ - ai^ f{asr - ai)^ 
ai + (1 - f)a2 



\ogp 



(5) 
(6) 



To achieve desired multiplexing gain r at high SNR (p 
00) the network must achieve a rate 7Z = r log p. The network 
is in outage if, min{/^^,/^^} < rlogp. For a given r and 
schedule / we can define the outage region 0(r, /) over 
channel realizations a = (ai, a2, <^sr). 



0{rj) = {a\ Y^l^<r} 



\ogp 



(7) 



The outage probability Vout is, 



Vout= / fd{<^l,C^2,(^sr) 

Jdeoirj) 

where /^i (o^i, 0^2, Q^sr) is the joint distribution of 
(ai,a2,Q^sr). 
Lemma 3.3: Proof see Appendix lAl 
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where < ai, 0^2 < 1, < agr < ^ and 
4 — 3ai — 2^2 



s{a) 



-{ 



T] H 

V- 



3ai - 
- 2ai 



Oi2 
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ai + ^2 < 1 
ai + 0^2 > 1 



(8) 



For a given listen-transmit schedule / the cut- set DMT 
bound is therefore given by, 



d(rj) = ^ min 5(a) 
deoirj) 



(9) 



To get the DMT upper bound we can optimize over all listen- 
transmit schedules, 



d(r) = min max 5(a) 

deoir) f 



(10) 



Note that this optimization is performed on a per realization 
basis, i.e. the optimal / depends on all channel realizations 
asrj(^i and a2. Therefore this corresponds to the global 
strategy discussed in Sec [III 

It is easy to see that the globally optimal schedule fgiob is 



one which sets 



Cs 



In 



fglob — 



a2 



{asr - <^l)+ + <^2 



This leads to the solution for d{r) given in ([T]). 

IV. ACHIEVABILITY: RELAYING SCHEME 

The "quantize-map" relaying scheme proposed in \T\ and 
|[8l| has been shown to be DMT optimal for the single an- 
tenna relay channel |5|. We show that "quantize-map" adapts 
naturally to the network described in Sec HI] and with some 
modification achieves the cut-set DMT bound. For the sake of 
completeness we include a short description of the scheme. 

A. Description of scheme 

S has a sequence of messages Wn ^ {1,2,..., 2^^}, n = 
1,2,... to be transmitted. At both source S and relay we 
create random Gaussian codebooks. S randomly maps each 
message to one of its Gaussian codewords and transmits it 
using T symbol times giving an overall transmission rate of 
IZ. Due to the half-duplex nature of the relay, it must operate 
using listen-transmit cycles. Relay listens to the first fT time 
symbols of each block i.e. X^. It quantizes Y^^ to Y^^ and 
then randomly maps it into a Gaussian codeword X^^ using a 
random mapping function /i?(Yi^). It transmits this codeword 
during the next (1 — /)T symbol times. Given the knowledge 
of all the encoding functions and signals received, D attempts 
to decode the message sent by S. 

B. DMT of Quantize-Map 

By Theorem 7.4.1 in (H, for any fixed listen-transmit 
schedule /, the quantize-map relaying scheme, uniformly over 
all channel realizations achieves a rate within a constant gap 
to the cut-set upper bound min{/^^,/^^} for that particular 
/. The random Gaussian code-book generated at source is 
independent of /. Also the code-book generated at relay 
depends on / only to determine the length of each codeword 
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Fig. 5. d{r) diocal{r) and duind{r) comparison for r > 1. 



The relay can generate a larger code-book with each code- 
word of length T. If the relay now chooses a listen-transmit 
schedule /, it can use the first (1 — f)T symbols of the 
codeword to compose X^^. The destination always knows the 
schedule / and hence can adapt its decoder accordingly. This 
construction allows us to claim that "quantize-map" achieves 
a rate within a constant gap of min{/^^,/^^} uniformly for 



each dynamic choice of / i.e. 



min{4^,4^} -n<R, 



quantize-map 



(hsr,K,K,pJ) (11) 



The constant in the above equation does not depend on the 
channel gains and SNR. At the order of DMT which assumes 
high SNR {p 00) the effect of n becomes negligible and 
hence we have the following theorem for achievability. 

Theorem 4.1: For dynamic listen-transmit schedules, the 
modified quantize-map relaying scheme as described above 
achieves the diversity multiplexing tradeoff of min{/^^ , /^^ }, 
where /^^ and /^^ are given by ([5])([6l). 

V. Achievability: Listen-Transmit Schedule 



In Sec IIII-AI the cut-set DMT upper bound was calculated 
for the globally optimal listen-transmit schedule fgioh- How- 
ever in a practical communication scenario global knowledge 
of instantaneous channel realizations may not be available 
at the relay. To account for this we defined the local and 
blind relaying strategies in Sec [III In this section we refine 
Theorem 13.11 to calculate DMT bounds for local and blind 
schedules. 

A. Blind Scheduling 

Theorem 5.1: For the low rate region i.e. r < 1, the blind 
scheduling strategy is DMT optimal. Additionally for 7^ > 2 
the blind strategy achieves the 2 x 2 MIMO DMT bound for 
r < 1. The optimal blind schedule for this region is fuind = | 



duindir) 



i{7/ + 2,4}-3r r<l 



(12) 



From ^ the DMT bound for blind scheduling can be 
written as, 

duind{r) = mdiX min s{d) (13) 

i.e. / is optimized without knowledge of channel realizations 
d. Solving this optimization for r < 1 yields Theorem 15.11 
This suggests that as long as cooperative multiplexing is not 
necessary i.e. desired rate 1Z = rlog(p) is such that r < 1, 
static scheduling at relay is sufficient to achieve the DMT 
upper bound, fuind = | turns out to be the optimal listen- 
transmit schedule for this region. 

For the high rate region (r > 1), the analytical solution for 
([T3]) is tedious to obtain. The optimization is convex and can be 
solved numerically. Fig. [5] shows a comparison between d{r) 
and duind{r) for r > 1. It can be seen that for cooperative 
multiplexing (r > 1) static scheduling is insufficient to achieve 
DMT upper bound. 

B. Local Scheduling 

Similarly, the DMT bound for local scheduling can be 
expressed as an optimization problem from (|9]). 



<^^oca^(^) = min max min s(d) 

(^sr f aua2eO{rJ) 



(14) 



fiocai can be optimized using knowledge of asr only. The 
DMT performance of local scheduling must be at-least as 
good as blind scheduling, therefore by Theorem lS.ll for r < 1 
diocai{r) = duind{r) = d{r). 

Numerical solution to ([T4l) for the high rate r > 1 region is 
shown in Fig. \5\ It can be seen that local scheduling performs 
better than blind, but for higher r] this advantage diminishes. 
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Appendix A 
Proof of Lemma[33] 

A. Marginal Distribution of agr 

fasri^) is calculated as. 



|2 < pC.,r-V] 



V[asr<a] = lim V[\hsrf < p''--''] 

p^oo 



B. Joint Distribution of ai and Oi2 

Note that g\ has a xl distribution, the marginal distribution 
of ai is given by, 

^^-2(l-ai)(-Q < < 1) 



Now for (0 < ai, Q;2 < 1) their joint CDF can be written as, 

Jx=0 



-2(l-x) 



x=0 



1 + 

Jx=0 

r p'V{\K\\r?<p-^^^-^) dy dx 
^y=o 



1) Case ai -\- a2 < I: 



Jx=0 Jy=0 



^ / ' 

Jx=Q 
^ ^3q!i+2q!2-4 

2) Case ai + 0^2 > 

Faua2{<^l^C^2) = / P' 
Jx=0 



' x=l — a2 

2Q!2+Q!1— 3 



f — 2Q!i+a2-3 

Jq!i,Q!2 r 

Since asr is independent of ai, 0^2 we get Lemma [33l 
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